Optical transition in self-assembled InAs/GaAs quantum lens under high hydrostatic 
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We present a simulation to characterize the dependence on hydrostatic pressure for the photolu- 
minescence spectra in self-assembled quantum dots with lens shape geometry. We have tested the 
physical effects of the band offset and electron-hole effective masses on the optical emission in dot 
lens. The model could be implemented to get qualitative information of the parameters involved in 
the quantum dot or the measured optical properties as function of pressure. 
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I. INTRODUCTION 

The study of the physical properties of nanostruc- 
tures is of great importance nowadays. Specifically, 
self-assembled quantum dots (QDs) grown by Stranski- 
Krastanov method^ have attracted great attention from 
the scientific community^'^. A lot of works have been 
devoted to the study of the electronic and optical prop- 
erties of QDs where their shape are modeled as quantum 
pyramids^i^iii^'SiiiSi^. It has been also considered QDs 
having lens shapei^i^^ii^^ii^ia^^. In this case, 
there is a rather extended opinion which states that the 
work done considering QDs with pyramids shapes can 
be straightforward applied to characterize the QDs with 
lens shape and therefore, the calculation made in pyra- 
mids QDs does not need to be re-calculated in quantum 
lenses. Nevertheless, different authors have shown the 
importance to consider the full lens geometry for a better 
description of the optical properties of the QDs'^^'^i^^iSi. 

Photoluminescence (PL) measurement under high hy- 
drostatic pressure has been proveed to be an effective 
tool for exploring the electronic structure and optical 
transitions in QDsS^^^^^^L^^i^^^^^^^L^. It is reported a 
significant decreasing in the hydrostatic pressure coeffi- 
cient (PC) for InAs/GaAs QDs in comparison with bulk 
material. Different experimental works have proposed 
some explanations: non- uniform Indium distribution in 
the QDs2S^ pressure dependence of the effective masses 
and confined potential^, and the necessity to include 
a nonlinear strain distributioit^^ and band-gap pressure 
dependence and finite barriers^!. From the theoretical 
point of view some works have been devoted to explain 
the obtained low value of hydrostatic PC. For example, 
elasticity theory in Ref. l25|. atomistic empirical pseudo- 
potential method in Ref. |28| and atomistic valence-force 
field method in Ref. [13 were implemented. In all of the 
former works, the QDs have been modeled as quantum 
pyramids. 

The motivation of the presente work is focused on the 



experimental observations of Ref. [2^[2^|30l[3ll . They 
have performed studies of low-temperature PL in InAs 
and InGaAs self-assembled quantum dots (SAQDs) em- 
bedded in GaAs matrix under hydrostatic pressure. 
From the pressure dependence of the main electron-hole 
transition energy, a pressure coefficient 20% smaller that 
the InAs case was reported. In Refs. [2^ and [29l some ar- 
guments could explain this difference in the pressure co- 
efficients: (i) the change of the band offset between InAs 
and GaAs, (ii) the strong tension in the InAs SAQD ma- 
terial and, (iii) the lowering of the electron and light-hole 
quantization energies due to the change of the effective 
masses with pressure. 

In Ref. [U was argued that the dots can be charac- 
terized by a lens symmetry with a base between 7.5 and 
10.6 nm of radius and ratio height/radius between 0.28 
and 0.4. We need to appoint that, in general, the shape 
and dimensions of the InAs/GaAs SAQDs are not always 
well determined and a degree of freedom remains. 

An attempt to describe the PL peak energies can be 
found in Ref. [s^ where the spatial confined geometry was 
modelled as a parabolic cylinder potential avoiding the 
specific geometry of the SAQDs. It is well known that dif- 
ferent geometries yield different results, as shown in Ref. 
[2^ when comparing the electronic properties between a 
cylindrical and a lens shape SAQD. In this paper, the QD 
will be modeled by a full lens symmetry with maximum 
height b and circular cross section of radius a with b < a. 
It is implemented a model which allows the inclusion of 
the finite barrier effect in SAQD with lens shape. The 
model allows to validate the importance of the change 
of the band offset and electron-hole effective masses as a 
function of pressure on the PL spectra. In our simula- 
tion the strain effect within the QD has been taking into 
account through an effective linear pressure coefficient at 
zero temperature which consider the built-in strain of the 
dot by assuming a simple two-dimensional strained layers 
model (see Table|T]and Refs. 32,35,36). The QDs consid- 
ered here are within the strong confinement regime and 
as first approximation the excitonic effects are ignored. 
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The remainder of the paper is organized as follows. 
Section [n] deals with the general trends of the simulation 
for including the hydrostatic pressure in QDs character- 
ized by a lens geometry with a short discussion of the 
model used to obtain the electronic spectrum. Section 
mil presents theoretical calculations for the system com- 
paring the effect of the hydrostatic pressure over different 
quantum lenses. A comparison with available PL data of 
the main transition energy as a function of P is presented 
in Sec. IIVI Finally, section |V] is devoted to the conclu- 
sions. 



II. BASIC RELATIONS 

The fractional change in volume in the InAs domain is 
given by 



SV/V = -3P{Sn + 2Si2), 



(1) 



where Sn = 1.946 x IQ-^ kbar'^ and = -6.855 x 
10^"' kbar~^ are the compliance constants^. The energy 
gap Eg at low temperature at the F point is proportional 
to the pressure according to 
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kT^ 



(2) 



where Eg"'' is the gap energy at T = and P = 0, a is 
the pressure coefficient and k and c are the temperature 
coefficients^. For the case of InAs/GaAs SAQD, the 
parameters used for both materials are shown in Table 
m With the inclusion of the finite stress it can be seen 
that the a coefficient is greater for the GaAs material 
than the InAs. Hence, the confined potential for the 
conduction and valence bands will increase with pressure. 
The variation of the effective masses was evaluated in 
terms of the fundamental gaps according to the Kane 
model and it was obtained^ 
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for light holes and 
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for the heavy holes. Here Eg is the gap energy given by 
Eq. ([5]), 71 and 72 are the Luttinger parameters, 5 is 
the spin-orbit splitting and nio is the free electron mass. 
The parameter Pq is the interband momentum matrix 
element between the conduction and valence bands. 

The PL lines in such QDs structures are proportional 
to the optical transition energies Ef4^^m^;Nh,mh{P) which 
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FIG. 1: Transversal section of a 3D lens Din in a semispherical 
medium Dout with contours C2 and Ci respectively. Contour 
Ci and C2 are separated a distance equal or greater than A 
along the perpendicular axis. 



are function of the pressure. The energy levels at a given 
band are characterized by the radial label N with ax- 
ial projection of the angular momentum with quantum 
number m. Then, in the framework of the parabolic band 
model, the transition energy EN^^rn^-Nh.rah means the en- 
ergy difference between the states {Nh: nih) and (iVe, rrie) 
in the valence and conduction bands, respectively. The 
Oscillator Strength (OS) is given by^: 



d{Ne,me;Nh,mh) 



^N^,m,:Nh.Tnh{r,r)dr 



(6) 



In the case of a SAQD with axial symmetry, it is obtained 
for the former expression d{Ne, me', Nh, mh) = for nif, ^ 
vrih and only the terms with me = rrih remains. The 
integration is over the whole domain D of the QD. 



A. Electronic structure 

Although previous theoretical studies in QDs with lens 
shape^i^i^ considered infinite potential wall, the finite 
value of the band offset is a fundamental parameter when 
considering the hydrostatic pressure, as can be seen from 
Eq. (21). Then, in our case, the problem for a finite bar- 
rier will be modelled including a lens-shape well potential 
with height Vo in a hard-walls semi-spherical region, as 
shown in Fig. [T] The semi-spherical region is divided in 
two regions. Din and Dout with confined potential zero 
and Vo, respectively, and the mismatch boundary condi- 
tions are considered at the interface. Assuming parabolic 
band dispersion for electron and holes, the solution of the 
effective mass equation for a lens geometry, with finite 
barrier, in an infinite surrounding medium can be ob- 
tained by minimizing the effect of the external boundary 
Ci over the wavefunction of the corresponding energy 
level under study. This can be achieved taken larger 
enough value of the distance A. The equation for the 
whole region D is: 



2 \m*{r,6,i 



V^] +V{r,e,(j))-^ ^ E-^, (7) 
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where 



and 
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The analytical solution of Eq. ([7]) is sought in the 
closed form 



i > 



(8) 



where the set of functions {^l""*} is a complete set of 
functions in the 3D domain D determined by the semi- 
sphere and its explicit representation can be found in 
Ref. [3^. With the former expansion the functions ^' 
satisfy the boundary condition of infinite barrier in the 
contour Ci because the set of functions does. On 

the other hand, Eq. ([7]) and the corresponding solution 
given by Eq. ([5]) are written for the whole domain D. It 
guarantees that the matching conditions at the contour 
C2 are also satisfied, but only at those points where it 
is well-defined the derivative of the wavefunction. This 
does not occur at the corners of the contour C2 where 
the differentiation with respect to the normal to the sur- 
face presents a discontinuity and, generally speaking, the 
problem is then not well-defined^. The obtained eigen- 
values constitute only an estimation of the solution for 
the real problem but we accept this solution as a bet- 
ter one than the values reported by the infinite barrier 
case. The present formalism has been applied in previous 
works^, but not explicit analysis has been done about 
the method used and the fulfillment of the matching con- 
ditions between the internal and the external domain. 

In order to visualize the results of the above de- 
scribed model, two different quantum lens geometries of 
InAs/GaAs have been considered. In Fig. [2] are shown 
the first five electronic levels as a function of the lens 
radius. Panel a) is for b/a = 0.91 while panel b) is for 
b/a — 0.51. The levels calculated within the hard wall 
model?2 are shown by dotted lines while the results ob- 
tained using the finite barrier model are described by 
solid lines. For the finite barrier case, it has been used a 
500 X 500 matrix and a value of A/R = 0.3 to guaran- 
tees that the absolute and relative errors introduced by 
the contour Ci and the cut-off of Eq. ([5]) in the values 
obtained for the energy levels are less than 10"'^. Ac- 
cording to the results of the Fig. [2l it can be seen that 
the infinite barrier is a good approximation for radius of 
the order of 20 nm or higher when b/a = 0.91. As b/a 
decreases, the values of the energy levels increase and the 
effects of the finite barrier become important. Then, it 
is necessary to consider the finite barrier effects to get 
better accuracy for the energy levels distribution at the 
same range of radius. 
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FIG. 2; First five electronic levels in the conduction band for 
an InAs/GaAs quantum lens as a function of the lens radius, 
a) h/a = 0.91, b) h/a = 0.51. The calculation is done taken 
A/R = 0.3 (sofid lines) for the finite barrier model. Dotted 
lines are the results considering the infinite potential barrier 
model^S. The conduction confined potential 14 is represented 
by dashed lines. 



TABLE I: Values of the parameters used in the calculation 
of the energy levels for the InAs/GaAs quantum lens under 
pressure. The symbol * means that the tension effect of the 
In As embedded in the GaAs matrix was taken into account. 





InAs 


GaAs 


El"' (eV) 


0.533''^* 


1.519" 


a (meV/kbar) 




10.8'= 


K (xlO-* eV/K) 


2.76" 


5.405" 


c(K) 


83" 


204" 


71 


20.4'' 


6.85'= 


72 


8.37'' 


2.1^= 


S (meV) 


380'' 


341'' 


2 (P^/mo) (eV) 


19'' 


12.9^ 



Ref. 44, '= Ref. 35, 
Ref. 45, " Ref. 46, Ref. 47 



" Ref. 43 
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III. RESULTS AND DISCUSSION 

Figure [3] shows the variation with pressure of the first 
energy levels in the conduction band for two lens con- 
figurations b/a = 0.91 [Fig. EJi)] and b/a = 0.51 [Fig. 
[5Jd)] with a — 9 nm. In both figures it is included the 
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FIG. 4; Oscillator Strength as a function of the pressure 
for electron- hole transitions {Ne,Nh;m) with Ne = Nh and 
m = 0. In solid line are shown the transitions from heavy 
holes states while in dashed lines from light hole states. The 
configuration lens is a = 9 nm with b/a — 0.91. In the inset 
the same radius and b/a = 0.51. 



FIG. 3: First energy levels in the conduction band as a func- 
tion of the pressure, a) b/a = 0.91 and a = 9 nm. b) 
b/a = 0.51 and a — 9 nm. In both panels the confined poten- 
tial energy it is shown by dashed lines. 



height of the potential barrier Vo — Vc a.s a. function 
of the pressure by dashed lines. The band offset of the 
strained InAs/GaAs quantum lens was taken, for the con- 
duction (valence) band, as 54% (46%) of the total band 
difference^. 

It is interesting to note the decreasing of the values of 
the conduction energies levels with the increasing of P, 
in spite of the increasing of the barrier height according 
to Eq. ^ and the decreasing of the dot radius following 
Eq. ^ . This effect is very robust for both configurations 
shown in Fig. [S] This fact is a direct consequence of the 
strong influence of the variation of the electron mass with 
P according to Eq. ([3]). The same will happen for the 
light hole (sec Eq. (g])). In the case of Fig. Ob), due to 
the strong lens deformation, only two electron states are 
present at P = [(1,0) and (1,1)] and two more appear 
at P w 20 kbar. 

In Fig. [5] it is shown the OS as a function of the pres- 
sure, for the first electron-hole states with m — 0. Each 
line is labelled according to the corresponding transition 
{Ne, Nh',m), which means a transition between the va- 
lence band states Nh to the conduction band states N^. 
Only those states with A^e = are shown, because they 
present the stronger values of the OS. The values of the 
OS corresponding to the transitions with Nh ^ N^ are 
very small and decrease with the pressure. There are 



shown by solid lines the transitions from heavy hole va- 
lence band and by dashed lines from light hole valence 
band, for comparison. The configuration used was a = 9 
nm and b/a — 0.91 while in the inset the same radius 
it is used, but a more confined region with b/a = 0.51. 
It can be seen that as the number A^ increases, the OS 
value decreases for a given pressure. Nevertheless, the 
OS for the light hole states have larger value than those 
corresponding to the heavy hole contributions. From the 
figure follows that as the ratio b/a decreases {b/a — 0.51 
in the inset) and keeping constant the lens radius, the 
value of the OS is reduced. This result is a direct con- 
sequence of the spatial confinement. Since the energy 
levels rise with the confinement (see Fig. [3]) the wave 
functions are less localized and the overlapping of the 
states with different effective masses decreases. Finally, 
it can be seen the increasing of the OS as P increases. 
Hence, we can conclude that at high pressure the main 
contribution to the observed PL lines comes from these 
states with A'e Nh- 

Taking in consideration the selection rules above ob- 
tained, in Fig. [5] is shown, in solid lines, the varia- 
tion with the pressure of the first allowed optical tran- 
sition energies for N^ — Nh and two InAs/GaAs quan- 
tum lens configuration. In Fig. [5] a) it is considered 
b/a = 0.91 while in Fig. [5] b) we have a stronger de- 
formed lens with b/a = 0.51. In both cases we used 
a = 9 nm, A/R = 0.3 and a 500 x 500 matrix in the 
diagonalization procedure. The dashed line is the varia- 
tion of the energy gap of the GaAs surrounding medium 
as a function of the pressure. In general, all transi- 
tions correspond to the electron-heavy hole states, ex- 
cept {Ne, Nih'jm) = (l,l;0)^j^, which corresponds to the 
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FIG. 5: First allowed optical transition energies of two 
InAs/GaAs quantum lens as a function of the pressure with 
a — 9 nm. a) b/a = 0.91. b) b/a = 0.51. Solid lines are ob- 
tained considering finite barrier while dotted lines represent 
the result with infinite barrier using an effective radius a^ff 
as a fitting parameter at P = 0. 
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FIG. 6: Ground state transition energy for a InAs/GaAs 
SAQD structure. Solid circles are data the reported in Refs. 
[2^ and [33. Solid line represents the theoretical calculation for 
a lens shape geometry with radius a = 6.4 nm at P = and 
b/a = 0.34. 



consideration that Fig. [5] b) corresponds to the strong 
confinement regime and the infinite barrier model break 
down. It the case of a quantum lens with higher values 
of b/a and radius (soft confinement regime), the energy 
levels are deep enough in the well. Then, the finite value 
of the band offset, and its dependence with the pressure, 
will have no appreciable influence on the transitions en- 
ergies and the model of infinite barrier will be suitable 
for characterize the system. 



IV. COMPARISON WITH EXPERIMENTAL 
RESULTS. 



electron-light hole state. As can be seen, the values of the 
transition energy for these states always increases with 
pressure. This behavior is due to the InAs energy gap de- 
pendence on P according to Eq. (l2|--'48_ ^j^g ^g^gg 
Fig. mb), the stronger confinement provokes an increases 
of the corresponding energy levels and at P = we are in 
presence of only two transitions, while transitions (1,2;2) 
and (2,2:0) appear for P > 20 kbar. 

On the other hand, each panel shows, in dotted line, 
the calculation using the infinite potential barrier model. 
Here, we introduced an effective radius fle// as a fitting 
parameter at P = for the ground transition energy 
(1, 1; 0). In Fig. [5] a) the fitting procedure is reasonable 
good for the electron-hole ground state (1,1;0) for all val- 
ues of P but it becomes worse for the excited ones. In the 
case of Fig. [5]b), the implementation an Ueff at P = 
does not work properly even for the (1,1;0) transition 
and it does not have sense for the (2,2;0) electron-hole 
state transition. This result can be understood taking in 



In spite of the limitation of the present simulation we 
explored its viability comparing with available experi- 
mental data. Figure [S] shows, in solid circles, the experi- 
mental transition energy for the electron-hole ground sate 
as a function of the pressure for an InAs/GaAs SAQD 
reported in Refs. [23 and [sO- We have calculated the 
ground state transition energy {Ng = 1, Nhh = 1; m = 0) 
in terms of the applied hydrostatic pressure. To fit the 
data in our calculation we fixed the ratio b/a — 0.34 inde- 
pendent of P and used the value of a = 6.4 nm at P = 0. 
The theoretical result is shown in solid line. Note that 
the value of b/a, but not the radius a, here em ploy ed 
match very well with the estimation given in Ref. l3ll 

Notice from the Fig. [5] that the value of the linear pres- 
sure coefficient is dE/dP ~ 8.3 meV/kbar. In connection 
with this, the variation of the band offset and the change 
of the effective masses with P are the responsible of the 
obtained value. Also, it is important to remark that the 
input used value of the linear pressure coefficient in Ta- 
ble [T] for the InAs is linked to the good match between 
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the simulation and the experimental data. If the typical 
bulk value of 100 meV/kbar for InAs is used the agree- 
ment goes down drastically. The consideration of the 
built-in strain of InAs/GaAs in the framework of two- 
dimensional layer model is one of key factor that allows 
the correct charactcirization of the fundamental optical 
emission as function of the hydrostatic pressure. 

V. CONCLUSIONS 

It was developed a model which allows us to include 
the pressure dependence of the volume, the band offset 
and the effective masses in the calculation of the energy 
levels at the F point of the Brillouin zone in SAQDs with 
lens shape geometry. It was outlined the importance of 
the involved parameters on the electronic states and on 
optical emission for SAQDs under hydrostatic pressure 
ranging between and 40 kbar. Also, we have delim- 
ited the validity of the infinite potential barrier model 
for the evaluation of the energy levels as a function of 
the hydrostatic pressure in quantum lens. 

The strong strain effect of InAs/GaAs on the gap en- 
ergy and pressure coefficient can be considered under 



two-dimensional layer approximation guarantying a reli- 
able description of the fundamental PL line as a function 
of the applied pressure. 

The present model is a qualitative one which needs to 
be improved in order to include the more complex effects 
present in such structures such as the real strain field in 
dots with lens symmetry, the peculiarity of the valence 
bands in III-V semiconductors and excitonic effects. 
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